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Abstract 

We propose a mechanism for the creation of cosmic string loops with dynamically sta- 
bilised windings in the internal space. Assuming a velocity correlations regime in the post- 
inflationary epoch, such windings are seen to arise naturally in string networks prior to loop 
formation. The angular momentum of the string in the compact space may then be sufficient 
to ensure that the windings remain stable after the loop chops off from the network, even 
if the internal manifold is simply connected. For concreteness we embed our model in the 
Klcbanov-Strasslcr geometry, which provides a natural mechanism for brane inflation, as 
well a being one of the best understood compactification schemes in type IIB string theory. 
We see that the interaction of angular momentum with the string tension causes the loop to 
oscillate between phases of expansion and contraction. This, in principle, should give rise to 
a distinct gravitational wave signature, the future detection of which could provide indirect 
evidence for the existence of extra dimensions. 



1 Introduction 



The existence of string loops with dynamically stabilised winding in a compact space was first 
demonstrated by Iglesias and Blanco-Pillado [I], using the Klebanov-Strassler geometry [2]. 
They considered strings at the tip of the throat, with geodesic wrappings in the S s which regu- 
larises the conifold singularity. Although they derived a lower bound for the angular momentum 
of a loop with a given number of windings - below which the windings became unstable - the 
result must remain of purely theoretical interest to cosmology as long as a specific mechanism 
for winding formation (and hence for the formation of the string angular momentum) is not 
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considered. The purpose of this paper is to propose such a mechanism, which leads to the for- 
mation of string configurations such as those investigated in [T] and to investigate the resulting 
string dynamics with specific reference to their cosmologically observable consequences. 

Assuming that a velocity correlations regime in the post-inflationary epoch as in [3] leads natu- 
rally to the formation of geodesic windings we show that the winding number (n), total energy 
{E) and angular momentum (I) of the string are specified precisely by the model parameters. 
That is, by the parameters which define the Klebanov-Strassler geometry (in this case specif- 
ically by the value of the warp factor ao) and those which determine the scale of the string 
network (a, ti). Substituting for I and n in the bound referred to above then demonstrates the 
stability of these windings, at least under the assumption that I remains approximately constant 
over small time scales after the moment of initial loop formation. By assuming also that the 
total energy of the string remains approximately constant (i.e. by neglecting the loss of E and 
/ via gravitational wave emission over small time scales), we then determine the equation of 
motion for the four dimensional string radius r(i), and solve to find a (generically) oscillating 
solution. Crucially we observe that the qualitative behaviour of the loop depends on the value of 
the warp factor ao with ag < 1/2 leading to an initial phase of expansion and Oq > 1/2 leading 
to an initial phase of contraction. The fixed point solution Oq = 1/2 is a static, non-oscillatory 
solution. In both oscillatory modes (initially contracting and initially expanding) we find that 
the period of the oscillation is inversely proportional to Oq, and proportional to the initial size 
of the loop {pdi). Following [3] we refer to these objects as non-topological cycloops|f|. 

The layout of the paper is then as follows: In Section 2 we briefly review the relevant Klebanov- 
Strassler background, focusing on the geometry of the conifold tip. In Section 3 we show how the 
assumption of a velocity correlations regime yields a dynamical model of winding formation after 
the end of inflation. Section 4 recaps the generic results of [1] which we then combine with the 
results of the previous section to give explicit expressions for the winding number, energy and 
angular momentum of the string. The equation of motion for the loop radius is then derived 
and solved in Section 5 and a brief summary of the main results together with a discussion 
of their cosmological implications and possible consequences for experimental observations is 
presented in section 6. Finally the two appendices at the end of this work deal with issues 
which arise within the text: Appendix I outlines the method of Eulerian substitution of the 
third kind which is used to integrate the differential equation involving r(t) and r(t) derived in 
Section 5. Appendix II gives a detailed description of the Hopf fibration of the 3-sphere, which 
is introduced briefly in section 2 and used throughout the following analysis. 

4 The term 'cycloop' was first coined in [2] to refer to cosmic string loops with smooth windings wrapping cycles 
in the internal space - as opposed, for example, to non-smooth, step-like windings which give rise to necklace 
configurations from a four dimensional perspective [4] [3 [6j . However in their original conception Avgoustidis and 
Shellard used it only to refer to windings which are topologically trapped. We therefore propose the term 'non- 
topological cycloops' to refer to string loops with dynamically stabilised smooth windings (in this case geodesies) 
around a simply connected compact manifold. 
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2 The Klebanov-Strassler Geometry 



The Klebanov-Strassler geometry is the canonical example of a background which resolves a 
conifold singularity in type IIB string theory. We refer the interested reader to the original 
paper [2] for a more thorough discussion. The essential point is that the conifold is the cone 
over an S 2 x S 3 base space. When we deform the conifold, the S 2 shrinks to zero size, and the ten 
dimensional metric factorises into the (warped) product of Mi^ x S 3 . In canonical coordinates 
the effective metric of the Klebanov-Strassler geometry at the tip of the warped throat takes 
the following form; 

ds 2 = c%q lu ,dx> i dx v + R 2 {dip 2 + sin 2 ^{d9 2 + sin 2 Odcf 2 )) (2.1) 



where r\uv is the usual four-dimensional Minkowski metric, a$ is the square of the warp factor 
such that < Oq < 1 and R is the radius of the three-sphere defined by, 

R 2 = bMg s a'. (2.2) 

Here M is the number of units of flux wrapping the internal space, which is also equal to the 
number of fractional Z)3-branes at the bottom of the throat, g s is the string coupling, l s = ya! 
is the fundamental length scale of the string and b is a numerical constant of order one. The 
relation between the size of the S 3 and the warp factor induced by the back reaction of the 
fluxes is, 

2-4/3 

where the constant e -4//3 is the deformation parameter of the conifold. Of course the relations 
(|2.2p and (|2.3p remain true in any coordinate system, and we will find it convenient to use 
the Hopf fibration of the three-sphere when considering the formation of geodesic windings. 
These are the kind of windings we expect to form in the presence of velocity correlations which 
impart an initial (constant) angular momentum density to each point along the string. In 
Hopf coordinates the S 3 is described as a one-parameter family of flat two-tori (to which it is 
topologically equivalent) [7] and the canonical metric (|2.ip reduces to a much simpler form, 

ds 2 = -cftq^dxPda? + R 2 (di? + d6 2 + cos## 2 ). (2.4) 

The Killing vectors also adopt a simple form, and are always parallel to the unit vectors in 
the ip, 9 and (^-directions. Fixing the value of the ^-coordinate such that 6 = 6q then selects 
a flat T 2 sub-manifold and windings which follow the Killing-directions in this manifold are 
necessarily geodesic in the full S 3 . However the choice of gauge in this respect is somewhat 
arbitrary, as we are free to choose 8q = 0. This simplifies both the resulting metric and the 
Killing vectors of the T 2 , the latter now being identical to the unit vectors in the remaining 
ip and ^-directions. Although it may be shown explicitly that the Lagrangian density C for 
a string loop with geodesic windings in the S 3 is cr-independent in any coordinate system, 
the simple form of the Killing vectors in Hopf coordinates allows us to more easily calculate 
L = j daC = 2tc£. A thorough treatment of the Hopf fibration of the three-sphere, and of the 
description of geodesic windings in both canonical and Hopf coordinates is given in Appendix 
II, along with a coordinate-independent geometric analysis. 
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3 A dynamical model of winding formation 



We now proceed to construct our dynamical model of winding formation: In the velocity corre- 
lations regime the number of windings per loop - in a loop of size r(U) = ati - is [3], 

n {ti) ~ — ^- (3.1) 

where ui is the fraction of the total string length which lies in the extra dimensions and is defined 
via [6], 

nR 

Ul ~ / -2 9 9^ T ( 3>2 ) 

V agH + n^it^ 

Substituting this back into the expression above gives a unique physical physical solution, 



n(U) ^ — • (3-3) 

We note that the condition ag < 1 ensures that n(tj) > OVij > 0. Importantly one notes that 
the number of windings increases linearly with time, with an overall coefficient modulated by 
the presence of the warp factor. 

Alternatively one can choose to solve for wi rather than the winding number; 



indicating that the magnitude of the warping imposes a physical constraint on the length of the 
string in the extra dimensions. We can also identify this quantity with the velocity of the string 
in the compact space. Imagine that the "end point" of the string at the horizon moves with a 
fixed velocity - < 1 in the extra dimensions. Then at time tj there will be approximately, 

n(U) ~ ^ (3.5) 

windings within the horizon (with c = 1) and the number of windings within a fraction a of the 
horizon is therefore, 

nft) ~ ^r- (3-6) 

This is equal to the number of windings per loop, for loops formed at time ti in the scaling 
regime. With this identification we see that 



Ul ~ v ~ y 1 - ag (3.7) 

and note the physical conditions < v < 1 and < uj[ < 1 are automatically satisfied by the 
condition < a§ < 1. Furthermore it is intuitively obvious that uji, v,n(ti) — > as ag — > 1 
because the ag = 1 solution of the Klebanov-Strassler model corresponds to Minkowski space 
in (9 + 1) dimensions. In this case there are no internal fluxes, and thus no compact extra 
dimensions implying that no windings can exist. 

5 Although this is not obvious from the formula (|2.3p . the rationale behind this statement is the following; if 
no fluxes exist to provide an effective potential with which to compactify the extra dimensions, there can be no 
back reaction on the ordinary four-dimensional Minkowski manifold. Hence Oq = 1 and R 2 — > oo automatically, 
leading to a flat six-dimensional space in place of the metric (|2.ip / (I2.4|l . Similarly the other dimensions of the 
bulk CY space are no longer flux-compactified. In such a scenario the formula (|2.3|l would not be valid as, by 
definition, it holds only for ajj < I . 
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We can also understand why ui, v — > 1 as <2q — > if we realise that windings are effectively 
correlations which can only form within the horizon. The horizon in the infinite dimensions 
advances according to the expression, 

d%(t) = a t (3.8) 

whereas in our background the horizon distance in Calabi-Yau space (S* 3 ) is simply, 

d c H Y (t) = t. (3.9) 

Although the expression Oq = is strictly unphysical (corresponding to an extremal horizon), 
the limit Oq — > corresponds to a situation in which the infinite dimensions are "closed off" so 
that the string exists only in the compactified space (hence uji = 1). In a time interval 5t, the 
value of limits the increase of the horizon distance in the infinite directions (but not in the 
compact space) via, 

<$d°° = a 5t (3.10) 

which places a limit on how fast the correlations can form. Strictly speaking, the horizon in the 
infinite direction advances by a distance ao5t when the horizon in the compact space advances 
by 5t. The end point of the string, which must of course move with resultant velocity v res = 1, 
must therefore cover a total distance in the compact space 5d given by, 



Sd = Jl- alSt (3.11) 

which limits the effective velocity of the string in the extra dimensions to, v ~ % = y/l — Oq. 
The parameter u>i is then given simply by the ratio of the string velocity in the compact space 
to the speed of light, oji ~ v ~ \J\ — Oq. Happily we find that the details of the compactification 
scheme determine v, coi and n(U) uniquely |j. 



4 Comparison with the results of Blanco-Pillado and Iglesias 

We wish to consider a string loop which has windings over the full S 3 . A previous study was 
initiated by pQ in which Euler variables were used to describe the S 13 as an SU(2) group manifold. 
This foliates the three sphere into a one-parameter set of flat tori. For fixed angle #o> the metric 
then reduces to that of a flat two-torus. The strings therefore only wrap a two-dimensional 
sub-manifold of the full three-sphere. We wish to generalize this result to consider windings 
over the full S 3 using the Hopf map. We use the following ansatz to describe a string loop, with 
general, non-specific windings around the S 3 

X M {a,t) = (t,r{t)sm{a),r(t)cos{a),z ,0,0,0,0,^(a,t),e(a,t),^(a,t)) (4.1) 

6 Note that in the preceding section it could be argued that, accounting for the effects of warping, the original 
formula for n(ti) in un- warped space (|3.1|l (which was taken directly from [3]) should be modified to give n(ti) ~ 
" l °°^° t '^ . However using this in conjunction with (|3.2[) gives n — u)i — as the only possible solution. As there 
are good physical grounds (outlined above) for believing that the identification cji ~ v ~ y/l — a§ (|3.7[) is valid 
we therefore choose to leave the formula H3.1|) unchanged, even in the presence of warped space. 
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where we have chosen our gauge so as to identify the world-sheet coordinate r with the time 
coordinate in the Lorentz frame of the loop, i.e. r — > t. In keeping with the physical scenario 
we are considering, we specify the ansatz completely so that i/j(a,t), 9(a,t) and 4>(a,t) describe 
geodesic windings. As previously stated, in Hopf coordinates the Killing vectors of the three- 
sphere take a particularly simple form and are parallel to the unit vectors (1,0,0), (0, 1,0) and 
(0,0, 1), so that geodesies windings are described by; 

ip(a,t) = In^a + tjj(t) 

9(a,t) = n e a + 9(t) (4.2) 
(p(a,t) = n^a + <p(t) 

where n^, rig, € 7L represent the number of physical windings present in each angular direction. 
Note that the factor of two in front of the term is simply an artifact of the coordinate system, 
resulting from the fact that the principle range of the V'-coordinate is twice that of 9 and <\> 0. 
Plugging ()4.ip and (|4.2p into the standard Nambu-Goto term of the D-string action (and ignoring 
the topological Chern-Simons term and other flux-dependent contributions) we have; 

S = -T x J dadt^a 2 {l - r 2 )(a 2 r 2 + R 2 s' 2 ) - a 2 r 2 R 2 s 2 (4.3) 

where s = ip + 9 + (ft and T\ is the tension of the .D-string. Rescaling the tension to absorb 
the string coupling allows us to also consider the F-string. In the warped deformed conifold, 
the tension spectrum of a general (p, q) string (in the large p or large q limit) is given by the 
well-known formula [U El QUI E] 

which we assume to be valid for p, q >> 1 and p « M in order for back-reaction effects to be 
neglected. 



In fact we may simplify our result even further by an appropriate gauge choice with respect to 
the angular coordinates. Since we know that geodesies of the compact space correspond to great 
circles on the S" 3 , we may set either of the winding numbers rig or to zero without loss of 
generality. The resulting string action with geodesic windings on the S 3 (in Hopf coordinates) 
is therefore, 

S = -Ti J dadt^/a 2 (l - r 2 )(a 2 r 2 + R 2 (2n^ + n^) 2 ) - a 2 r 2 R 2 (ij; + 0) 2 . (4.5) 

The resulting constants of motion are then; 

^ dL dL , , . 

E = —q-L l = —q' (4.6 
dq q 

where q G {t, a, ip, 9, 0}. The first expression is simply the Hamiltonian, and the second corre- 



7 The principle ranges of the angular coordinates are 0<^<27r, O<6 , <7r,O<0<7r. 
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sponds to the total angular momentum of the string in the compact directions |f| Again using 
our ansatz (|4,ip - (|4,2|) we see that these expressions become: 

E 2nT a 2 (a 2 r 2 + R 2 (2n^ + n e ) 2 ) 



ajftl - r 2 ){a 2 r 2 + R 2 (2n^ + n^f) - a 2 r 2 R 2 (ijj + </>) 2 

a 2 r 2 R 2 (2n i) + n e )(ip + e) 
a 2 (l - r 2 )(a 2 r 2 + R 2 (2n^ + r^) 2 ) - a 2 r 2 R 2 (iP + 4>) 2 

Now Iglesias and Blanco-Pillado [I] have shown that, for a loop which is stationary in 3 + 1 
dimensions (i.e r = 0), the energy of the string configuration is minimised precisely for, 

r 2 = (4-8) 

and 

s 2 = (^ + 0) 2 = § . (4.9) 

These results are obtained by first re- writing s 2 = (ip+6) 2 in terms of I, r, R and s' 2 = (2n^+ng) 2 
and then substituting into the expression for E (so that E = E(r, I, R, s')), finally minimising 
with respect to r. Technically this gives the first condition (|4.8p . and the second condition 
is then obtained by further substitution into the expression for I. 



Now any dynamical model we construct for the formation of geodesic windings and for the 
motion of the string after loop formation must be consistent with these general results. At first 
sight we see that our model suggests v 2 (ti) = s 2 (ti)R 2 ~ (1 — Oq), which does not correspond (in 
general) to the energy minimisation condition v 2 = s 2 R 2 = Oq. In fact these two conditions only 
coincide for the specific value = 1/2, where the velocity in both compact and non-compact 
dimensions is v ~ 1/V2- We would not expect such a string configuration to undergo time 
evolution under the influence of its own internal dynamics, although it may still "shrink" via 
the loss of mass-energy (and angular momentum) due to gravitational wave emission and this 
possibility is discussed in section 6. We therefore conclude that if the value of the warp factor 
is exactly ao = l/y/2, the energy of the string configuration will be automatically minimised 
from the moment of loop formation i.e. the initial radius of the string loop r(ti) and the initial 
angular momentum l(ti) will be related via r 2 = Z/27tTiciq, and the velocity of the string in the 
compact space will be v ~ a/1 — = ao = l/\/2. 

However when a^ > 1/2, the string velocity in the compact space is too small to provide enough 
angular momentum to "match" the radius of the loop i.e. the angular momentum required 
for a cycloop of radius r to minimise the energy of its configuration. Alternatively one can 

8 Note that although the string "rotates" around the S 3 , no centripetal force is acting upon it. The internal 
(compact) dimensions are parameterised in terms of the angular variables ip, 8 and (f>, and so the motion through 
the S 3 may be measured in rad x [t] -1 . As the effective radius of "rotation" for any point along the string is 
simply the radius of the three-sphere, multiplication by R converts this "angular velocity" into the "true velocity" 
of the string. However, even here we must be careful - as the string has no internal structure, the "velocity" of 
the string parallel to itself (in this case parallel to the geodesic windings) is not clearly defined. It is therefore 
possible (in principle) for v(t) ~ s(t)R > 1 though this does not violate causality due to the boost invariance of 
the string along its length. 
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understand this as a string loop with fixed energy, changing configuration in order to minimise 
the surface to energy ratio. Hence from the arguments in Iglesias, we expect that ag > 1/2 
implies I < 27rTiagr 2 . Similarly if ag < 1/2 this should imply that the initial angular momentum 
of the loop exceeds the optimum value for a loop of that size and we find the converse result, 
namely I > 2 , KT\a\r 2 . In a full string theory compactification, the warp factor is exponentially 
small and therefore this would appear to be the dominant string channel. However we will 
consider a more phenomenological approach in this paper, and consider a range of values for 
the warping. Furthermore, as stated above, these results should hold true for any physically 
viable dynamical model. It is therefore worth testing the theory developed in Section 2 to ensure 
consistency in this matter. In short Iglesias' and Blanco-Pillado's arguments [T] regarding the 
energy minimisation condition (for all r, I) ought to be consistent with our own dynamical model 
of l(ti) and r(ti) outlined above. 



Because we are dealing with geodesic windings, we may always re-define our coordinate system 
so as to identify the variable n from (|3. 1 j) with the variable s' , both of which represent the total 
number of physical windings in the compact space - hence n ~ 2n^ + ng. We now plug in the 
expressions for nfc) (gTB, v 2 = s 2 R 2 = ag (JO) and r(U) = aU into P~T|) (with r(U) = 0) to 
find the resulting terms □ 

E( ti ) = 2VTT, ao{aU) 



V^o 2 

l(U) = 2irT ia 2 (ati) 2 . (4.10) 



Hence we see that the second part of the energy minimisation condition (|4.9p implies the first, 
and vice- versa, as expected for consistency. However, as noted above, in general we have v 2 (ti) = 
s 2 (ti)R = (1 — Oq) from ()3.7p which is not equal to ag unless ag = (1 — ag) = 1/2 giving the 
constant values 

E = E(ti) = 2-irT^ati) 
I = l(ti) = 7rT l (at t ) 2 . (4.11) 

Under such special circumstances we would not expect the string configuration to evolve in time 
due to its own internal dynamics, though the emission of gravitational waves due to accelerated 
motion of the string, and the resulting shrinkage of the loop radius r must still be accounted for 
as mentioned previously. 



Considering the more general case (a 2 + 1/2) and substituting n{t { ) (HD}, v 2 (U) = s 2 (ti)R 2 = 
(1 — ag) (|3.7|) and r(tj) = ati into (|4.7|) - keeping f(U) = - gives us 

E(U) = 27rTi(ati) 
l(U) = 2ttTi(1 - a 2 ) (ati) 2 (4.12) 

9 Note that we may could have substituted n(ti) = a^(odi) 2 / R 2 in place of the usual expression (|3.1|l . taking 
advantage of the fact that v 2 ~ a^R 2 in this case. This leads to the expression E(ti) = 2nT\ x 2a1,(ati) which is 

,2 2 



equivalent to (|4.10|) for a — (1 — a ) = 1/2. We therefore see that the total energy is split half and half between 
the rest mass of the loop in warped Minkowski space, and the kinetic energy due to motion in the S 3 . 
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which implies that 

l(U) > 2irT ia 2 {ati) 2 a 2 < - 

l(U) < 2irT 1 a 2 (at i ) 2 a 2 > -. (4.13) 

We find that the total energy of a cycloop with radius r(ij) = ati is independent of <Zq. At 
first glance this seems nonsensical: the value of the warp factor determines the velocity in the 
compact space at the moment of loop formation via v(ti) ~ s(ti)R ~ yl — Oq, which in turn 
determines the initial number of loops via n(ti)R ~ v(ti)r(ti) ~ y/l — a^ati). A cycloop moving 
with greater velocity in the compact space would therefore have a greater number of windings 
than a slower moving string with the same radius in the non-compact directions. Consequently 
an increase in the kinetic energy of the loop would seem to go hand in hand with an increase in 
the total rest mass. However although this is clearly true in un-warped space, we must remember 
that the presence of the warp factor also reduces the four dimensional energy density via the 
effective tension T\ = c^T\. Equation (|4.12p suggests that even though a smaller warp factor 
implies a greater rate of winding formation and a greater kinetic energy for the windings in the 
compact space, the would-be increase in the total energy of the cycloop is completely off-set by 
the reduction in four dimensional energy density. 

The question then remains: what happens if the energy minimisation conditions are not au- 
tomatically satisfied at the moment of loop formation ij? This is of course equivalent to the 
question: What happens dynamically when either Oq < 1/2 or > 1/2? Intuitively we would 
expect that if l(ti) > 2i:Tia^r{ti) 2 (oq < 1/2), the radius of the loop would rapidly expand 
introducing a non-zero r(t) term for t > ti lj. Physically this corresponds to the conversion 
of kinetic energy from the motion of the string in the compact space, into rest-mass energy in 
four dimensions. The conservation of angular momentum also suggests that any fraction of lifi) 
'lost' in this process (whatever this proportion may be) is carried away by the gravitational 
radiation produced by the the expanding loop, though as a first approximation we may neglect 
this. Hence we must allow for the most general case by including the explicitly time-dependent 
term r = f(t) in the expressions for E and I (|4.7p . whose derivative r(t) we expect to be initially 
positive for an expanding loop (i.e. r(ij) > 0). Needless to say, we must also re-introduce a 
time-dependent velocity term v(t) = s(t)R for t > ti, whose derivative v(t) we expect initially 
to be negative in this case (v(U) = s(ti)R < 0). 

Similarly if l(ti) < 2irTia^r{ti) 2 (oq > 1/2) we expect the opposite process to occur - with 
rest-mass energy of the loop being converted into kinetic energy in the extra dimensions. For 
t > ti we again introduce the extra dynamical terms r(t), whose derivative r(t) we now expect 
to be initially negative (r(U) < 0) and v(t) = s(t)R whose derivative v(t) we now expect to be 
initially positive (v(ti) = s(ti)R > 0). Again we face the possibility that a significant proportion 
of the initial angular momentum of the loop at the moment of formation will eventually be 
lost through the emission of gravitational radiation during dynamical evolution. As a first 
approximation however, we will assume the loss of angular momentum via gravitational wave 
emission to be negligible, taking I « l(U)Vt > ti. We will also assume that the total energy 
lost via gravitational wave emission during the dynamical evolution of the loop is negligible i.e. 

10 Note that this inequality is strict. At t = U we still have that r(ti) = 0. If this were not the case then the 
energy minimisation condition (|4.8|) would itself be different. 
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that E ~ E(ti) Vt > tj. What drives the evolution in this case is not changes in the energy of 
the system to "match" the conditions, but changes in the conditions to match the given energy: 
that is, the mutual and inter-dependent evolution of r(t) and s(t) toward a loop configuration 
which meets the energy minimisation criteria (|4.8p - (|4.9|) . 

The approach outlined above has the added advantage that in both cases we may assume that 
the number of windings remains fixed. As we shall see in the next section, the stability of the 
extra-dimensional windings places a lower bound on the value of /. If I remains constant, all 
that is required to ensure stability of the windings throughout dynamical evolution towards 
the minimum energy state is that the string have sufficient angular momentum to stabilise its 
windings at the moment of loop creation. In the next section we will demonstrate (following 
the analysis in PQ) that the stabilisation of windings places a bound on I, with more windings 
requiring a larger angular momentum to remain stable. Thus we see that, if the value of I were to 
change significantly during the evolution process, the dynamics of the loop may be considerably 
more complicated with windings "falling off" the S s as the loop expands/contracts. 



5 Loop dynamics after formation 



We now investigate the stability requirements for the extra-dimensional windings (as mentioned 
above) both generically and in light of our specific dynamical model. Happily we find that the 
rate of winding production in the model we have put forward ensures the stability of the ansatz 
(|4.ip - (|4.2p at the moment of loop creation for all possible formation times ij. We then leave this 
result on one side and consider the loop dynamics in each of the two different regimes (l(ti) > 
2-KTia^r{ti) 2 , Oq < 1/2 and Z(tj) < 2-KTia^r{ti) 2 , Oq > 1/2) discussed above on the assumption 
that the loops retain their windings during the evolution towards an energy-minimising state 0, 
which is equivalent to the assumption that I = l(ti) for all t > 

By introducing a small perturbation in one of the bulk-space directions perpendicular to the 
S 3 it is possible to show, from the resulting expansion for the ten-dimensional action, that the 
string configuration (|4.ip - (|4.2p is stable [1] (this corresponds to moving the entire string "up" 
from the very bottom of the throat by a small amount). The ansatz (|4.2p (with ng = 0) also 
implicitly assumes that the motion of the string in the compact space is parallel to the direction 
of the windings, so it is not physically meaningful to 'perturb' the string in either the ip or in- 
directions. We may however investigate the effect of perturbing the string in the ^-direction in 
order to determine the stability of the winding configuration. Turning on a small perturbation 
58 then results in the following perturbation of the Lagrangian pQ 

a 2 

SL = (l/2TrT ia 2 + R 2 {2n^ + n e ) 2 ) 58 2 - {l/2^T ia 2 - R 2 {2n^ + n e f) 59 2 (5.1) 

11 Note that we will also find, in the next section, that the loop does not remain stable at the energy-mininimising 
configuration. However it is still true to say that it evolves from the initial radius towards such a configuration. As 
we will show, the loops actually 'overshoots' its own energy-minimising configuration due to the non-zero velocity 
of the radial coordinate at that point, leading to an oscillating solution. 
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which results in the stability condition 



I > 2-KT 1 alR 2 {2n il + n e f (5.2) 

In other words, the total angular momentum must satisfy this bound (note also the strictness 
of the inequality) in order for the number of windings {2n^ + rig) to remain stable. Again 
identifying n ~ 2n^ + rig and substituting for n{ti) (j3.3|) and l(t{) (|4.10p . we may investigate 
the stability of the windings in our dynamical model at the moment of loop formation. The 
resulting condition simply reduces to, 

a 2 < 1 (5.3) 

which is automatically satisfied by the definition of al in the Klebanov-Strassler geometry. Thus 
we see that the stability condition is satisfied and that all windings are stable at the time of 
loop formation for all ti and for all physical values of a,Q, R 2 and a. 



We now introduce a non-zero time dependent term r(t) for t > ti in (|4.7p which initially satisfies, 

r(U) > al < - 

r(U) < a 2 >^ (5.4) 
and a time-dependent velocity term v(t) = s(t)R which initially satisfies, 

s(ti)R = v(U) < al < - 

s{ti)R = v{U) > al > - (5.5) 
and which also satisfy the following boundary conditions, in the limit t — > ti, 

s{ti)R = v(U) ~ ^/l-og. HU) = 0. (5.6) 

From now on then we may assume that the number of windings remains constant from the 
moment of loop formation, n = n{ti) Vi > ti, and attempt to determine the corresponding 
dynamical evolution of the loop in the Minkowski directions. Using I = l(ti) and E = E(ti)Vt > 
h we then have 



l = l(ti) = 2ttTi(1 - al)(ati 



271-71(1^1 - a'l(ati)r 2 (t)s(t)R 



^a 2 (l - r2(t))(ogr 2 (t) + (1 - a^)(at l ) 2 ) - a 2 R^(t)s 2 (t) 

and 

E = E{U) = 2 7 rTi(at i ) 

27rT ia 2 (a 2 r 2 (t) + (1 - aDiaU) 2 ) 



(5.7) 



(5.8) 



Vag(l - f2(t))(ay(t) + (1 - aDiaUY) - a 2 R^{t)s^{t) 
Re-arranging (|5.7|) then gives, 

alr 2 s 2 R 2 = ag(l - ajj)(l - r 2 ){ati) 2 (5.9) 



11 



and re-arranging (|5.8p gives, 

2 2 . 2r 2 _ a 2 (a 2 r 2 (t) + (1 - ag)(c^) 2 )[(l - r 2 )(atj) 2 - a 2 Q (a 2 r 2 (t) + (1 - a§)(afr) 2 )] 
a r s R — ^ i 5 - 10 ) 

so that equating the two expressions above yields the following non-linear first order differential 
equation in r(t) 

r 2 + j^r 2 + (-1 + 2a 2 (l - a 2 )) + (1 - a 2 f{aU) 2 ^ = 0. (5.11) 

We note that the constant terms involving <2q and the terms involving powers of r form a perfect 
square, so that this equation may be re-written as, 

j*- 1+ (A- r+ ( 1 -<®W ) 2 = Q (5.12) 
\{ati) r J 

It is then explicitly clear that there exist two critical values of r at which r = i.e. at which 
the expansion of the loop (at least momentarily) comes to a halt. These are, 

r c x = (ati) (5.13) 

and 

(1 - # 



r C 2 - 2 



^(aU) (5.14) 



Although it is not possible to show this directly without the explicit form of the solution r(t), 
we know that the first of these values must correspond to the boundary condition r(ti) = 
which we imposed when calculating E(ti) and Z(tj), as well as when determining the energy- 
minimisation conditions (|4.8|) - (|4.9|) . However the second of these values is intriguing as it does 
not correspond to the minimum energy condition. We can tell immediately therefore that the 
dynamical evolution of the loop will not lead to a steady energy-minimising state, and we may 
may instead expect a solution which oscillates between the two values r c % = (ati) (the initial 

radius of the loop) and r& = <z°\ aii). 

a 



To see if such a solution is consistent with the physical arguments above i.e that a value of 
<2q < 1/2 leads to an initially expanding loop, Oq > 1/2 leads to an initially contracting loop 
and that Oq = 1/2 leads to a static loop of radius r(t) = (ati)\/t > t{, we must now ask 
two questions. Firstly, which is greater, the initial radius or the second critical value? This will 
determine whether or not the loop initially expands or contracts. And secondly, which is greater, 
the second critical value or the radius corresponding to the energy-minimisation condition? This 
will determine whether the behaviour of the loop is in accordance with our assumptions. 

Whether r c \ = r(ti) < (>)r C 2 therefore depends on whether (1 — ao)/ a o > or equivalently 
whether (1 — Og) > (<)«§• This in turn depends on whether ag < (>)l/2 with ag < 1/2 

12 Alternatively of course we may substitute either of the expressions (|5.10p or H5.9|) into the original string 
action (|4.5p and then determine the Euler-Lagrange equations. The resulting equations must necessarily have the 
same solution as (|5.11|) but the method adopted here which utilises the string constants of motion is far simpler. 
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implying that the loop must expand from its initial value r(ti) = (atA to r = — %°\ ati) > (atA 
and a 2 , > 1/2 implying that the loop must contract from its initial value r(ti) = (ati) to 
r = ^(0<«). 



The second question may be answered as follows: By equating our expression for l(ti) (|4,12p 
with the second part of the energy- minimisation conditions (|4.9p . we see that the critical radius 
corresponding to the fulfillment of this condition r m i n may be written in terms of a 2 , and the 
initial radius r(ti) = (ati) such that, 



v/r ^(at.) (5.15) 



on 



Whether r c2 < (>) r min then depends on whether (1 — Oq)/(2q < (>) — — — - which is equivalent 
to the condition (1 — afy/a^ < an< ^ reduces to a 2 , > (<)l/2. 

We therefore see that a 2 , < 1/2 implies that the loop expands from its initial value r(tj) = (ati) 

/l_ a 2 

towards the radius corresponding to the energy- minimising configuration r m j n = a(1 ° (o^) > 
r(ij) but overshoots it and continues expanding to the second critical value r C 2 = ^ %°\ ati) > 
r min- Similarly a 2 , > 1/2 implies that the loop contracts from its initial value r{t{) = (ati) 

fl—a? 

towards the radius corresponding to the energy-minimising configuration r m j n = 

r(ti) but overshoots it and continues contracting to the second critical value r C 2 = — (ati) < 

a 

r m i n . In both cases, of course, we expect the loop to oscillate back and forth between r c \ = r(t{) 
and r C 2 but our findings are consistent with the assumptions about the initial behaviour of 
the loop expressed in the boundary conditions (|5.4|) - (|5.6|) . We also note that for Oq = 1/2, 
r(ti) = r C 2 = r m i n = (ati) and the loop remains static in four dimensions. 

Having determined the general behaviour of our solution in each of the three cases Oq <=> 1/2, 
we now return to equation (|5.1ip in order to find the explicit form of r(t). By making the 
substitution y = r 2 (|5.11|) may be re-arranged to give, 

dt = ±\—= dV (5.16) 
2 V -ay 2 + by - c 

where we have defined the coefficients a, b and c to be explicitly positive for Oq in the range 
< Oq < 1 so that, 



(«^) 2 

6 = 1- 2ag(l - ajj) 

c = (l-alf(atif (5.17) 

We also note for future reference that the discriminant of the quadratic A = b 2 — Aac (given 
below) is positive for all < a 2 , < 1 except Oq = 1/2 for which A = 0. 

A = b 2 - Aac = 1 - 4o§(l - ag) = (1 - 2a 2 ,) 2 (5.18) 
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Integrating (|5,16p then yields the following explicit expression for t in terms of y, where K is 
simply the usual integration constant and A, B are the real roots of the equation —ay 2 +by—c = 
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t 



1 

tan 

\ a 



' B-y 
y-A 



which may be re-arranged to give y = y{t), 

B + At&n 2 (y/a(t + K)) 



y(t) 



1 + tan 2 (Va(i + K)) 



Now the roots of the quadratic —ay 2 + by — c = are given by, 



y 
y 



(i-4) 2 



■{ati 



(5.19) 



(5.20) 



(5.21) 



but choosing A = {ati) 2 and B = — {ati) 2 does not allow us to fix K in order to satisfy the 



boundary conditions. We must therefore choose B = {ati) 2 and A 
(|5.2Up before setting y{U) = {ati) 2 to determine, 

K = -U. 

Our final solution for r{t) is then given by taking r = +y/y giving; 



{ati) 2 in equation 



(5.22) 



r{t) 



{ati) 2 + ±^{aU) 2 tan 2 (J^{t - U)) 



which may be re- written in terms of basic trigonometric functions. 



(5.23) 



r{t) 



c2 



1 - 2a 2 . \ 

u 1 sur 



hit-U) 



{at 



(5.24) 



It is clear that for a 2 , 7^ 1/2, the solution (|5.23p / (|5.24p oscillates between r{ti) = {ati) at 
t = mr^^- for n E {0,1,2,...} and r, 



'^{aU) at t = (2m + l)§^, m £ {0,1,2,...}. 

We also see that for a 2 , = 1/2 we have A = B (A = 0) and r{t) = {ati) Vi > ti as expected. It 
is also tedious but straightforward to twice differentiate (|5.23p in order to verify the boundary 
conditions ()5.4[) . The relevant quantities are, 



r{t) 
r{t) 



;i - 2a 2 ) 



sin(x) cos(x) < 1 + 



1 



+ 



1 - 2af 



{ati 

cos 2 {x) — sin 2 (x) 



M sin 2 (:c) 



(5.25) 



) sin 2 {x) cos 2 (x) \ 1 + 







1 - 2af 



0, 







^ J sin 2 (2:) 



13 This solution was obtained by performing a Euelerian transformation of the third kind. A brief sketch of this 
method is given in Appendix I. 
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where we have defined x = j^hj(t — U), from which it can clearly be seen that as t — > t j we have 

r(ti) = and rfo) = ^7jS* , and hence f(t») > for < 1/2 and rfc) < for > 1/2 in 
accordance with the boundary conditions (|5.4p . 



Figure 1 below illustrates the qualitatively different behaviour of the solution (|5.24p for different 
values of Oq. The three curves show the three different types of dynamical evolution that a loop 
(formed with initial radius r(ti) = ati) may undergo in the ag < 1/2, ag = 1/2 and Oq > 1/2 
cases respectively. It is also worthwhile plotting the behaviour of v(t) ~ s(t)R for the same test 
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Figure 1: This figure illustrates the behaviour of the solution (|5.24|) in the three qualitatively 
different regimes. For convenience we have chosen ti = 1, a = 0.5 for all three curves and set 
Oq = 0.4, 0.5 and 0.6 respectively. 



values of Oq, a and U which is shown in Figure 2. The formula v(t) may be obtained by taking 
either (|5.9p or ()5.10|) and substituting for r(t) from (|5.24p and r(t) from (|5.25p . though we will 
omit the explicit expression here for the sake of brevity. As expected, the behaviour of v(t) (or 
rather v(t)) clearly satisfies the boundary condition constraints (|5.5|) . Even though, for any given 
model, the values of Oq and a are constant; string loops are continuously formed throughout 
the history of the universe, meaning that ij will be a continuously changing parameter. In each 
of the oscillating regimes we would therefore expect to see a spectrum of oscillation periods, 
with smaller loops formed at earlier epochs oscillating more rapidly between their initial and 
maximum radii than larger loops formed at late times. This is illustrated for the o,q > 1/2 
regime in Figure 3. The corresponding behaviour of v(t) ~ s(t)R is shown in Figure 4 using the 
same test values for the parameters. 
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Figure 2: Figure illustrating the behaviour of v(t) ~ s(t)R in the three qualitatively different 
regimes. Again we have chosen ti = 1, a = 0.5 for all three curves and the values = 0.4, 0.5 
and 0.6 respectively. 
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Figure 3: This figure illustrates the behaviour of r{t) for loops formed at three different epochs 
in the Oq > 1/2 regime. For the sake of convenience we have fixed clq = 0.4 and a = 0.5 for all 
three curves and set ti = 1, tj = 1.2, and tj = 1.5 respectively. 
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Figure 4: This figure illustrates the behaviour of v(t) ~ s(t)R for loops formed at three different 
epochs in the ag > 1/2 regime. As in Figure 3 we choose to set Oq = 0.4 and a = 0.5 for all 
three curves and consider tj = 1, ti = 1.2, and t, = 1.5. 

6 Discussion 

We have argued that a velocity correlations regime in the post-inflationary epoch leads naturally 
to the formation of string loops with geodesic windings in the compact space. For strings at 
the tip of the conifold throat of the Klebanov-Strassler geometry we were able to show that the 
quantities which determine the dynamical evolution of the loop (i.e. the initial winding number 
n(ti), energy E(ti) and string velocity /angular momentum in the compact space, s(ti)/ l(ti)) 
are uniquely determined by the parameters ajj, a and t^. Crucially these windings were found 
to have sufficient angular momentum in the compact directions to remain stable, even after the 
string chops off from the network to form a loop. 

The interaction between the tension and the angular momentum in the compact space was found 
to play a significant role in the dynamical evolution of the string, including - perhaps surprisingly 
- the evolution in four dimensions. By assuming energy loss (and angular momentum) via grav- 
itational radiation to be negligible, we determined equations of motion for the four-dimensional 
radius r(t) and the string velocity v(t) ~ s(t)R, which we believe to be valid over small time 
scales after the moment of loop formation, t = ti. 

We found that the qualitative behaviour of the string depends crucially on the square of the 
warp factor, < Oq < 1, with ag < 1/2 leading to an oscillatory solution characterised by 
an initial phase of expansion, whilst a§ > 1/2 leads to oscillations with an initial contracting 
phase. In each case the string was seen to oscillate between it's initial radius r(ti) = (ati) and 
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a secondary critical value denned by r C 2 = - — ^(aiA with the period of oscillation given by 

a o 

rp _ (ottj ) 
a 

As noted above, in the two oscillatory regimes it is the interaction of the angular momentum 
with the string tension which "drives" the dynamical evolution converting kinetic energy into 
rest mass during expansion (with the inverse process occurring in the contracting phase). The 
string is seen to evolve towards a static, minimum energy configuration where I = Ina^Ti and 
s 2 = Oq/-R 2 , but is unable to satisfy these two conditions simultaneously at a point where 
r(t) = 0. By contrast, in the cJq = 1/2 scenario we find that the energy minimising conditions 
are satisfied simultaneously at the moment of loop formation. In this case the tension exactly 
offsets the effect of angular momentum so that the string remains static at it's original radius 
r(ti) = (ati) and is non-oscillatory. 

The meaning of the term "small" above is somewhat ambiguous, but it seems reasonable to 
assume that our solution will provide a valid approximation over at least one full oscillatory 
cycle of the loop; that is, over a time period At ~ T = For time periods At >> T our 

initial analysis must be extended to include the effects of gravitational wave emission on E(t) 
and l(t) (or equivalently on r{t) and s{t)). We may expect the qualitative effect of energy and 
angular momentum loss to be relatively simple, as long as the string retains sufficient angular 
momentum for the extra-dimensional windings to remain stable. In this case it is likely that the 
loss of E and I due to gravitational wave emission will simply act to damp the oscillations of 
r{t) and s(t). What is unclear however is whether or not the damping coefficient itself is likely 
to be time-dependent; For example, it is possible that smaller oscillations lead to greater rates 
of emission per unit length (as the rate of acceleration r{t) is higher in this case) so that the 
damping itself increases with time (for an individual loop). 

In the case that the windings eventually become unstable however (as l{t) drops below the 
threshold for ensuring their stability) the string dynamics are likely to become complicated, and 
it is not clear whether the process of winding contraction (i.e. of windings "falling off" the S 3 ) 
may even be accommodated within an analysis which uses an ansatz of the form (|4.ip to describe 
the string configuration. This is because the coordinates r(t)sina, r(t)cosa, vp(o~,t), 9(a,t) and 
</>(<7, t) are treated as independent variables with respect to the determination of the equations 
of motion. We are therefore unable to take account of the continuously connected nature of the 
string when string sections "move" from one direction to another (e.g. in "falling off" the S 3 to 
form part of the four dimensional rest mass of the loop) . 

The present analysis could of course still be improved by accounting for the effects of gravita- 
tional wave emission under the assumption that that loops retain their windings, which would 
indeed be valid up to the point where l(t) drops below the critical value (which may also be 
calculated). Such an improved analysis could proceed as follows: one could compute the stress 
energy tensor for an oscillating loop and look for solutions with this as a source to the Einstein 
equations. This should allow us to estimate the rate of loss of E and I via gravitational wave 
emission and, as stated, we expect this to produce a damping term in the equations for r(t) and 
s(t). 
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Calculating the emission spectrum for an oscillating loop would also be of immense practical 
interest. The gravitational wave signature of such a loop, whose self-oscillation is caused by 
the presence of angular momentum in the compact space, may differ significantly from that 
of a string whose oscillations, though superficially similar, do not result from self-interaction. 
Gravitational wave emission from loops oscillating with period to ~ L _1 (where L is the loop 
size) have been intensively studied in four dimensions |16t [T71 [TBI [T9] . However in such cases the 
loop is not undergoing genuine phases of expansion and contraction, but rather experiencing 
"wiggles" of a size comparable to it's own length. Although Weinberg [20] has has shown that - 
in an FRW universe - the power of a weak, isolated, periodic source (to lowest order in G) may 
be given by a single formula regardless of the exact nature of the source, it is not immediately 
clear that this should hold in extra-dimensional scenarios. Additionally such sources (i.e. loops) 
have no angular momentum to shed in the process of emission. In fact even if we were to 
study loops whose self-oscillation was due to their "rotational" motion in Minkowski space (c.f. 
|21| 122]). the angular momentum which would be carried away by gravitational radiation would 
be very different from that lost via emission from oscillating cycloops 0. The detection of such 
a signature by the future generation gravitational wave detectors could then provide indirect 
evidence of the presence of compact dimensions. We hope therefore to be able to provide such 
an analysis in the near future. 

Although the full string theory compactification favours exponentially small values of <Zq, it is 
also interesting to note that the possibility of obtaining evidence from the gravitational wave 
signature of wound strings exists even in the case of static loops (i.e. the cJq = 1/2 case). 
These loops, which form automatically with a configuration which meets the energy minimising 
conditions, look just like ordinary string loops from a four-dimensional perspective. However 
they contain an " unseen" angular momentum which is not directly manifested in their dynamical 
evolution in Minkowski space. As stated above, we expect that even "standard" four dimensional 
string loops may undergo periodic oscillations with uj ~ creating ripples in space-time and 
giving rise to gravitational waves. Such fluctuations along the length of the string - though 
not in the total four dimensional length itself - would likely still occur in this case (though it 
is possible that the existence of the angular momentum term may lend a certain 'rigidity' to 
the circular string configuration, making it more resistant to deformation) but the string must 
also now shed it's angular momentum. The gravitational wave signature of even a non-self- 
oscillating loop in the extra-dimensional scenario is therefore also likely to differ significantly 
from the standard case of an un-wound string, though perhaps to a less significant degree than in 
the self-oscillating case. However the possibility of the indirect detection of compact dimensions 
from cosmic strings therefore remains, even in the absence of self-oscillating loops. We hope also 
to be able to offer an analysis of this interesting possibility in a future letter. 

It remains for us to outline some of the possible limitations of the analysis we have attempted: We 
have assumed throughout the present work that the velocity correlations regime leads naturally 
to a) geodesic windings and b) movement of the string parallel to itself (i.e. along the geodesies). 
However these assumptions may be questioned. The rationale for adopting such an approach 
(which simplified the resulting analysis considerably) was simple - we expect velocity correlations 
to impart a constant angular momentum density to each point along the string. Therefore we 

14 We note also that one possible further extension of the current analysis would be to consider a wound string 
with rotational motion in both the compact and Minkowski directions. 
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expect each "point" (or infinitesimal string segment) to travel along a geodesic curve, both 
before and after loop formation. However it has not been proved that each point along the 
string traversing a separate geodesic, leads to a winding configuration that is itself geodesic. 
Likewise it does not necessarily follow that the resulting motion of the string as a whole is 
parallel to itself. In general, we may expect that the geometry of the internal space plays a role 
in determining the exact nature of the winding configuration and the resulting string motion. 
Moreover we have not included the contribution from the Ramond-Ramond (RR) sector, which 
in principle will couple to the string. This charge term is ultimately crucial for distinguishing 
between cosmic strings and cosmic super-strings. 

Whilst, in principle, motion of the string perpendicular to its length may easily be accounted 
for (see end of Appendix II) the most significant problem in the analysis of string loops with 
non-geodesic windings arises from the resulting u-dependence of the expressions for E and I. 
We believe that the present analysis may therefore be improved by a more thorough investi- 
gation of the winding process itself, and though we expect our expressions for n(ti) (|3.3p and 
v(ti) ~ s(ti)R (|3.7p to remain valid (in the KS case), the appropriate string ansatz may well be 
more complicated. Although it seems unlikely that the qualitative behaviour of the string will 
differ significantly than described in the scenarios above, finding analogous (and quantitatively 
different) results may be extremely difficult. But ultimately exact quantitative predictions will 
be needed for any comparison with future experimental data, and we believe such a project to 
be worthwhile. 
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Appendix I: Eulerian substitution of the third kind 

As noted previously, the quadratic equation in the integral defined by (|5.16|) and (|5.17|) has 
a discriminant A = b 2 — Aac which is everywhere non-negative for values of Oq in the range 
< Oq < 1. We may therefore evaluate the integral using a Eulerian transformation of the third 
kind, which takes advantage of the fact that there exist two real roots in y. We proceed in the 
following way: Let A and B be the two real roots of the quadratic equation —ay 2 + by — c = 0. 
We then define a dummy variable u such that, 
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V ' —ay 2 + by — c = (y — A)u 



(A-l) 



In general of course we have that, 



ay 2 + by- c = a(y - A)(B - y) 



(A-2) 



so that combining (|A-ip and (|A-2|) implies, 



a{B -y) = (y- A)u 2 



(A-3) 
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or equivalently, 

aB + Au 2 .. A . 

y = — — — A-4 

a + u A 

and 



4= = (A-5) 

Using the second solution for y as function of u and differentiating the expression above then 
gives, 

2du f (aB + Au 2 \ A ,. n , 

dy = — -j \ - — — )+A\u (A-6) 

Finally, substituting (|A-ip . (|A-4|) and ()A-6p into the right hand side of (|5.16p gives 



2 J ^-ay 1 + by- c 

±1 f d y 



± 
2 



2 7 (y-A)« 

1 / a+u 



2du ( ( aB+Au 2 \ , /i 1 

tilt 



a + n 2 



(A-7) 



Using the standard integral f ^if^i = tan 1 (|) and (|A-5[> then gives equation ()5.19p . 



Appendix II: The Hopf fibration of the three-sphere 



As stated in the introduction, we chose to use the Hopf fibration of the S 3 when considering 
geodesic windings, as this allows the metric and Killing vectors to be written in a particularly 
simple form viz f|2.4[) / (|4.2j) . Choosing windings which wrap only Killing directions in the S 3 then 
leads to manifest a independence in the constants of motion (|4.7p . The ansatz corresponding 
specifically to geodesic windings is the linear ansatz (|4.2p , and the choice of coordinates naturally 
reflects the SO (3) symmetry of the internal space. 



We will now discuss the origin of the Hopf fibration from a geometric point of view. However it 
will also be useful to discuss the origin of the canonical coordinate system from a similar per- 
spective. By comparing the two (equivalent) descriptions, we hope to make clear the advantage 
of using the former coordinate system in the present work, but also to demonstrate explicitly 
the coordinate independence of our results. 

One very natural choice of coordinates on the <S 3 is the so-called 'canonical parameterisation'. 
In this case the group element of the S 3 manifold is given by, 

g _ e (i/2)[x<n +2/0-2 +za 3 ] ^ (B-l) 
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where cxj, i G {1, 2, 3} are the usual Pauli matrices. Writing the usual Cartesian coords x, y and 
z in terms of polars r, 9, eft the line element then becomes, 

ds 2 = -^Trddgg- 1 } 2 ) = dr 2 + sin 2 r{d6 2 + sin 2 Qd<\> 2 ) (B-2) 

which is exactly that given in f)2. 1 1) if we identify ip = r and multiply the whole metric by i? 2 . 

By contrast, and following Iglesias and Blanco-Pillado [I], we chose to parameterise the S 3 using 
Eulerian variables. In this case we regard the S 3 as an SU(2) group manifold with element g 
given by, 

g(if>, 9, <f>) = e iW/2)<Tl e i{9/2)cT2 e iW2)cT3 (B-3) 

where Oi are the Pauli matrices as before. The group invariant metric (i.e. the metric on S 3 in 
these coordinates) can then be written in the following form 

ds 2 = -hrr{[dgg- l f) = ^ + cos##) 2 + ^(d9) 2 . (B-4) 

It may then be seen that, in these coordinates, the sub- manifolds where 9 = 9q is constant 
correspond to flat two-tori with metrics given by, 

ds 2 = ^(dip + cos O #) 2 (B-5) 

This shows that the Hopf fibration corresponds to describing the three-sphere as a one-parameter 
family of flat two-tori (to which it is topologically equivalent). Again we follow Iglesias [1] in 
choosing #o = in order to simplify the metric as much as possible. It is of course also necessary 
to re-scale the metric so that 1/2 — > R 2 as we are dealing with a physical S 3 of radius R. 



Let us now consider an arbitrary string embedding which is a function of both space and time, 

X l (a, t) where i £ {0, 1, . . . 10}. The general - coordinate independent - expression for the string 
Lagrangian is then, 

L = -T x al j da^j (1 - r 2 )(r 2 + a^ 2 R 2 W) - a^ 2 R 2 r 2 P + a^R^Q (B-6) 

where 

W = X>,i) 9li pT(M))X'>,t) (B-7) 
P = X l (a,t)g^(X(a,t))Xi(a,t) 

Q = (XV^^M))*'^*)) 2 

- {X\a, t) 9lj (X(a, t))X' J (a, t))(X' l (a, t) 9ij (X(a, t))X^(a, t)) 

where a dash indicates differentiation with respect to a and a dot represents differentiation with 
respect to t as usual. In the static case the term W(a) = X' l (a)gij(X(a))X ,:J (a) corresponds 
simply to the length of the wrapped string on the S 3 . Demanding that = one simply 
obtains the geodesic equation, 

X"\a) + T) k X'i{a)X' k {o) = (B-8) 
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where Y l - k are the usual Christoffel symbols. Thus minimal (i.e. geodesic) windings on a general 
manifold (in our case this is just S 3 ) with metric gij is enough to guarantee a independence of 
the winding terms. 

The same statement still holds true if we allow time-dependence into the string ansatz, X l (a, t), 
as long as for each instant in time the latter still satisfies the geodesic equation. Physically this 
means that the wrapped loop only ever evolves along geodesic curves. A simple example on the 
S 2 would be a great circle passing through the north and south poles and rotating about an axis 
through the poles. Thus we require our embedding function to satisfy, 

X' H (a,t) + T) k X IJ (a,t)X' k (a,t) = (B-9) 

Now, with time-dependence we also have a non-zero kinetic energy term of the wrapped string 
inside the square root factor in our Lagrangian, 

P(a,t) = X'faQgijiXfatyXifat). (B-10) 

Demanding now that ^=0we find, 

X H (a,t) + Y) k X j (a,t)X' k (a,t) = (B-ll) 

Geometrically this equation states that the velocity vector, X 1 , is preserved under parallel 
transport along a geodesic curve that is wrapped by the string. One can equally interpret this 
equation as saying that under parallel transport of the tangent vector to the geodesic curve, 
X' 1 , along a curve whose tangent vector is X 1 , the former is preserved. Mathematically this is 
simply the statement that V^' = Vyy = 0. 

The nice consequence of demanding the above is that that the winding term is not only a 
independent (if we take geodesic wrapping) but also time independent, i.e. ^5- = 0. Thus a 
independence of the kinetic function P guarantees time independence of W . This is an important 
statement. 

Let us now consider what the situation above implies for the case where the wrapping is over 
flat sub-manifolds as in the case of Iglesias, where the T l - k all vanish. The geodesic equation 
(|B-9p and (|B-lip trivially imply X" 1 = 0, X H = 0, and therefore the solution of these equations 
is, 

X l (a,t) =X t + n i a + u i (t). (B-12) 

This provides the origin of the linear a anasatz of Iglesias. Now apart from the winding and 
kinetic terms mentioned above there are additionally the two R 4 terms appearing inside the 
square root in the Lagrangian in Q(a,t). The a independence of the second term follows from 
our previous results, because this term is simply WP. The a independence of the first term 
above follows after a simple calculation making use of (|B-9p and (|B-11|) . In the case of wrapping 
along a flat sub-manifold one sees that Q vanishes identically. 

However what is important is that a independence of all the relevant terms in the Lagrangian 
is guaranteed simply by requiring the strings to wrap geodesies and that the velocity vector is 
preserved under parallel transport along this geodesic. Thus we see that L = J daC = 2ir£ as 
stated in the introduction. 
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We now have the possibility of working in any coordinate system, though the choice of explicit 
string wrapping ansatz is constrained by the requirement of solving ()B-9|) . (|B-1 1|) if we require 
geodesic windings. In general explicit solutions to the geodesic equations are hard to come by 
but for very symmetric spaces - like spheres - they are known Ej . 



However let us now compare the form of these explicit solutions, in canonical coordinates, to 
the simple linear ansatz in (|B-12p . In order to determine the appropriate ansatz for geodesic 
windings in canonical coordinates we must first calculate the Killing vectors in this coordinate 
system. We know that the Killing fields corresponding to the SU(2) rotations of the S 3 generate 
isometries of the above metric - therefore, if we think of the S 3 as being embedded in R 4 where 
X, Y, Z, U represent the four Cartesian coordinates, 

X 2 + Y 2 + U 2 + Z 2 = 1, (B-13) 

the three independent isometries are generated by 

J\ = Jxu + Jyz, Ji = Jxz + Jyu, -h = Jxy + Jzu (B-14) 

where Jzu = Zdjj — Udz is the generator of rotations in the Z — U plane etc. The above Ji 
clearly generate an SU (2) algebra. Then the three Killing vectors are, 

k\ = (-U,-Z,X,Y); 4 = (-Z,U,-Y ) X); fc* = (-Y, X, Z, -U) (B-15) 

which, in fact, define an orthonormal basis for the SU (2) Lie algebra. However unlike the Killing 
vector fields of the two-torus, those of the whole S 3 are in general coordinate dependent. 



Since we work with coords ip, 0, <f> rather than the embedding coords X, Y, U, Z one can re-express 
the above Killing vectors in terms of k^ijj, 6, 4>) with i = 1 ... 3 of the 5 3 . To obtain these just 
consider the left action of rigid group elements 

9a = d eaa (B-16) 

on the group element g of SU (2) written in terms of canonical coordinates described above. 
Explicitly then we have that, 

gn = cos(V0 + i sm(ip) cos(#) (B-17) 

gi2 = ism(ip)sin(9)e~ l ' t ' 

g 2 i = i sin(V>) sm(8)e l ' t ' 

922 = cos(V0 — i sin(^) cos(#) 

Expanding out to linear order in e in the gi and the reading of the infinitesimal variations 
Sip, 59, 5cj> and by equating 

5X i = ekid j X i (B-18) 



lj For general wrappings of the string around the full S however, the ansatz will certainly never be linear in a 
and nor will the a and time dependence factorize in an additive way as it did in the flat space case, no matter 
what choice of coordinates we make. 
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we can then read off the components of the killing vector k l a for each of the isometries induced 
by left action with group element g a , a = 1, 2, 3 of SU (2). This gives, 

h = \sm(9) cos(», sin(» - cot(^) cos(0) cosU), - cot ((9) cos((p) - cot(^) (B-19) 

sin(t') 

/c2 = [— sin(#) sin(<^>), cos((/>) — cotf-0) cos(6>) sin(0), — cot(#) sin(<^>) — cot(tp)— — -^-1 

sm(t/) 

fc 3 = [cos(0),-cot(^)sin(0),-l] 

One can check (its a bit tedious!) that the above killing vectors are orthonormal with respect 
to to the canonical metric on S 3 

KmA = S ab (B-20) 

Although, using the above results, it is possible therefore to show the a independence of the 
Lagrangian for geodesic windings, even when working in canonical coordinates, the resulting 
expressions are incredibly complicated. If we now return to our coordinate independent de- 
scription - armed with the knowledge that the Lagrangian density must be a independent for 
geodesic windings - we find that we can in fact reproduce most of our results without reference 
to a specific coordinate system. 



Looking at the effective potential (setting all time derivatives to zero) we find, 

V = aln^r 2 + clq 2 R 2 W (B-21) 

and we now know that the second term is just the warped length squared of the wrapped string 
on the S 3 . The tangent vector to a geodesic x H (where a dash here refers to differentiation with 
respect to the affine parameter along the curve) can always be defined to have unit length ie 
x ll gijX l:/ = 1. In our case the string wraps closed curve geodesies with different winding numbers 
in general. A nice example of closed curve geodesies on the <S 3 are the integral curves whose 
tangent vectors are the three Killing vectors fe*,a = 1,2,3 that generate the SU(2) isometry 
group of S 3 (see above). Thus it is natural to take the following ansatz for our static wrapped 
string 

*g(*) = E = E n aK(X(a)) (B-22) 

a a 

where we have used the fact that the tangent vectors to the SU(2) generated geodesies are just 
the Killing vectors k\ (and where the subscript G implies that the embedding ansatz satisfies 
(|B-9|) and (|B-lip ). The Killing vectors not only have unit length with respect to the canonical 
metric on the S 3 but they are in fact orthonormal. Using this it is easy to see that W = n a n a , 
which is a constant and so is consistent with our previous analysis. 



The above makes sense because it is known that flows of the Killing vector fields on S 3 induced 
by rigid SU (2) rotations are actually geodesies. So our wrapping ansatz would be guaranteed to 
be 'minimal' in the same sense that the linear ansatz of Iglesias for winding around the torus - 
because again the flows of the Killing vectors in that case are geodesic circles around the torus. 

Notice that in equation (|B-22p we have written the expression for X i (a) and not X l (a). The 
latter can be obtained in principle by integration. In the simple case where the Killing vectors 
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are coordinate independent (e.g. Abelian isometries, shift isometries), integration directly gives 
us Xq = Xq + an a k l a which we found in the case of flat sub-manifolds discussed earlier. The 
SU(2) Killing vectors however are not constant, so integration is non-trivial and will lead to 
non-linear dependence on a in general. 

It is a general result that if a Killing vector field has constant length, then it will generate 
geodesies along the manifold. In our case we know that the k l a (X) have unit length with respect 
to the canonical metric on the S* 3 . Therefore we can be sure that equation (|B-9j) will be solved 
by X l% {a) ~ k l a {X) on the three sphere according to the ansatz in (|B-22p . One can see this more 
clearly by taking the normalisation constraint gijk l a k 3 h = 5 a b and differentiating this with respect 
to some arbitrary vector field Y. The resultant Killing equation gives us V^fcj, = 0, which more 
generally implies 

Vn ak n b k b = (B-23) 
which is just the restatement that n a k l a is a geodesic. 

Since we are interested in the dynamics of such an embedding we now need to extend our 
solution to incorporate time dependence. This means that we need to preserve the geodesies 
under time evolution. To ensure this, let us use the SU(2) transformations, since these are 
isometries preserving the canonical metric, and therefore we map geodesies to geodesies. The 
modified ansatz for the embedding functions should now read 

X*(<r,t) = ^n a ki(X(cr,t)) 

a 

X\a,t) = ^2x a kl(X(a,t)) (B-24) 

a 

where we have introduced new variables A a £ R. The string modes will therefore be wrapped 
along any curve defined by n a k l a , and this curve will then evolve in time via the second geodesic 
equation (with tangent vector X a k l a )- This means that the functions in the action will take a 
simplified form 

W ->• ^nana (B-25) 

a 

P -> J^A a A a 

which are both constant, and A a is related to the average speed of the string along the S 3 . 

Now it n a and A a are such that the X H and X % are parallel, then the dangerous R 4 terms will 
vanish from the action. If these vectors are not parallel then there is no cancellation, which 
means that the wrapped geodesic must have a perpendicular component along the winding 
direction satisfying X l X' J gij = 0. 

How does this affect the resultant action constructed in equation ()B-6j) ? We can define a unit 
vector pointing parallel to the winding direction in the usual manner. This therefore allows us 
to split the velocity into components X % = X^ + X l j_ . After some manipulation we find that the 
Lagrangian simplifies to give 

L = -Ta 2 ^(l -P- a^^XigijX^ir 2 + a~ 2 R 2 W) - a- 2 R 2 r 2 X^ gij X^ (B-26) 
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where the i? 4 terms have canceled as already advertised. Physically this makes sense, since the 
perpendicular modes lie in the normal bundle and therefore contribute to the transverse boost of 
the string, much like the velocity in the Minkowski directions. The net effect is an enhancement 
of the relativistic 'gamma' factor. 

In the simplest case where we neglect the transverse modes, we see that the resultant energy 
and momentum become 

a%(r 2 + a^ 2 R 2 W) 

E = 2vrTi 



(1 _ /i2)( r 2 + a^R 2 W) - dQ 2 R 2 r 2 P 

R2r2XaHa , (B-27) 

'(1 - r 2 )(r 2 + a^ 2 R 2 W) - a^ 2 R 2 r 2 P 

where we have defined I in the same manner as Iglesias 

I = ^X'\ (B-28) 

For the case where there is no velocity in the Minkowski direction, we can write the energy as 
a function of I and minimise it to obtain 

r - ~ (B - 29) 

which is the generalisation of the r ~ \fl dependence obtained in Iglesias. Moreover the velocity 
at the minimal radius is given by 

X\ = XaXa = Jr- (B-30) 

Generally we should consider the case where there is non-zero X l ± , we will see that solutions 
that minimise the energy whilst having non-zero I require us to set f = X l ± = - which is easily 
understood since non-zero velocity in these directions only ever increases the energy through 
enhancement of the 'gamma' factor. 

In the general case we see that the above expressions can be written in the modified form 

(r 2 + aQ 2 R 2 n a n a ) 



E = 2ttTi 

^(1 _ r 2 _ afR 2 X i ± g l3 X{)(r 2 + afR 2 n b n b ) - afR 2 r 2 \ b \ b 

I = 2itT x ; tf^Kng = (B31) 

1 _ r 2 _ a fR 2 X i ± g lJ X{){r 2 + afR 2 n b n b ) - a Q 2 R 2 r 2 \ b \ b 



References 



[1] J. J. Blanco-Pillado and A. Iglesias, JHEP 0508, 040 (2005) [arXiv:hep-th/0504068] . 
[2] I. R. Klebanov and M. J. Strassler, JHEP 0008, 052 (2000) [arXiv:hep-th/0007191] . 



27 



[3] A. Avgoustidis and E. P. S. Shellard, JHEP 0508, 092 (2005) [arXiv:hep-ph/0504 049|. 



[4] T. Matsuda, arXiv:hep-ph/0601014 



[5] T. Matsuda, JHEP 0604, 017 (2006) [arXiv:hep-ph/0509062] . 
[6] M. Lake, S. Thomas and J. Ward, arXiv:0906.3~695l [hep-ph]. 

[7] I. Bengtsson and K. Zyczkowski, Geometry of Quantum States (2006), CUP, ISBN: 978-0- 
521-81451-5 (hardback) 

[8] J. H. Schwarz, Phys Lett. B 360, 13 (1995) [Erratum-ibid. B 364, 252 (1995)] 
[arXiv:hep-th/9508143] . 

[9] H. Firouzjahi, L. Leblond and S. H. Henry Tye, JHEP 0605, 047 (2006) 
|arXiv:hep-th/0603i6T] . 

[10] S. Thomas and J. Ward, JHEP 0612, 057 (2006) [arXiv:hep-th/0605099| . 

[11] H. Firouzjahi, JHEP 0612, 031 (2006) [arXiv:hep-th/0610130] . 

[12] T. Matsuda, JHEP 0505, 015 (2005) [arXiv:hep-ph/0412290] . 

[13] T. Matsuda, JCAP 0604, 005 (2006) [arXiv:hep-ph/0509064] . 

[14] A. Avgoustidis and E. P. S. Shellard, Phys. Rev. D 71, 123513 (2005) 
[arXiv:hep-ph/0410349] . 

[15] V. Berezinsky and A. Vilenkin, Phys. Rev. Lett. 79, 5202 (1997) jarXiv:astro-ph/9704257l . 

[16] T. Vaschpati, and A. Vilenkin, Phys. Rev. D30, 2036 (1985) 

[17] C.J. Burden, Phys. Lett. 164B, 277 (1985) 

[18] D. Garfinkle and T. Vaschpati, Phys. Rev. D36, 2229 (1987) 

[19] B. Allen and E.P.S. Shellard, Phys. Rev. D45, 1898 (1992) 

[20] W. Weinberg Gravitation and Cosmology (1972), Wiley, New York, ISBN: 978-81-265-1755- 

8 

[21] R. Durrer, Nucl. Phys. B328, 238 (1989) 

[22] A. Vilenkin and E.P.S. Shellard, Cosmic strings and other topological defects (2000), CUP, 
ISBN: 0-521-65476-9 (paperback) 



28 



